Hopping-disorder-induced effects upon the two-Magnon Raman Scattering 

spectrum in an Antiferromagnet 
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Two-magnon Raman scattering intensity is obtained for two-dimensional spin- 1/2 an- 
tiferromagnets in presence of hopping disorder. A consistent mode assignment scheme is 
prescribed which essentially establishes a correspondence between the hopping-disordered 
and the pure system. It is seen that a minute amount of disorder in hopping leads to a good 
agreement of the Raman lineshape experimentally obtained for the copper-oxide insulators 
such as La2Cu04. Also it is observed that a considerable asymmetry with respect to the two- 
magnon peak appears and the spectral intensity persisting much beyond 4 J (where the joint 
magnon density of states peaks). This has been argued earlier by us [S. Basu and A. Singh, 
Phys. Rev. B, 53, 6406 (1996)] to be arising due to highly asymmetric magnon-energy 
renormalization because of a cooperative effect arising from local correlation in hopping 
disorder. 
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Recently it was put forward by Nori et al. jj] that the exchange disorder caused by zero-point lattice vibra- 
tion can account for various experimentally observed features of the two-magnon Raman intensity lineshape. 
In their approach the exchange disorder is looked upon as an interaction between the lattice excitation (i.e. 
phonons) and the spin excitations (i.e. magnons) which induces changes SJij , in the exchange integral J. An 

estimate of the mean square deviation ^(— j 2 -) 2 ^ is made and to make a comparison with the experimental 

data the corresponding value is considered to be as high as 0.5 (in units of J). It has been suspected Q that 
such a large variation in J appears physically unreasonable and questioned that whether more moderate val- 
ues of (—f 2 -) are sufficient to achieve agreement with experiments. Essentially what one obtains is that the 



coupling of phonons to the spin excitations leads to the broadening of a two-magnon peak, the asymmetry of 
the lineshape about its maximum towards the higher energy regime and appearance of scattering intensity 
in forbidden geometries like Ai gi B 2g etc. 

Earlier, we have investigated, [[| the effect of lattice fluctuation induced hopping disorder in the context 
of Mott-Hubbard AF where the hopping term include random terms. The Hubbard Hamiltonian in this 
case can be written as, 



H 



E 

<U>." 



(t + 5tij)(a^a jiCT + h.c.) + U 



nitnu 



(1) 



zexp 



2(7 



(2) 



where the random terms dUj's for each bond are chosen independently from a gaussian distribution and 
the distribution width y/a measures the strength of disorder. There we have examined both perturbatively 
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and using the exact eigenstate method - effects of hopping disorder on magnon energies, their wavefunctions 
and on the density of states (DOS). In both schemes, the spin- wave energies are obtained at the Random 
Phase Approximation (RPA) level. We have obtained a strong renormalization in energy of the high energy 
magnon modes, arising from locally correlated hopping, which results in appreciable one-magnon DOS well 
beyond the maximum spin-wave energy 2 J for the pure system. This result has got important implications 
in explaining the experimentally observed Raman spectrum. 

To obtain the Raman scattering intensity we need to consider the imaginary part of the full two-magnon 
propagator which in turn requires the knowledge of spin-wave energies. For a finite size system these 
energies are obtained numerically at the RPA level in presence of disorder using the exact eigenstate method 
discussed in detail elsewhere Q . Following a prescription of the mode assignments discussed in section III, 
it is possible to obtain the symmetry factor required for the B\ g scattering geometry (in Appendix A, we 
discuss the calculation of the symmetry factor in the B\ g geometry corresponding to the pure case). It is 
observed that even for a sufficiently small strength of disorder, a = 0.01, (much smaller than the values 
considered by Nori et al. though the Raman spectrum on the lower energy side is weakly affected, 
there is a substantial change in the higher frequency regime and a good agreement is obtained with the 
experimental data. This confirms our earlier statement H that if indeed the randomness in hopping due 
to zero-point lattice fluctuation is playing a key-role in explaining the long standing puzzling feature of the 
Raman lineshape regarding the asymmetry, then it essentially arises from the highly asymmetric magnon- 
energy renormalization. 

II. CONFIGURATION AVERAGING OF THE TWO-MAGNON PROPAGATOR 

The Fleury-Loudon (FL) Hamiltonian which represents the interaction between the photon and the spin 
pairs can be written as, 

H R = A^(E inc i).(E sc i)S(r).S(r + 5) (3) 

r.S 

Where the E.;„ c and E sc are the incident and the scattered electric field vectors, S is the unit vector connecting 
nearest neighbour (NN) sites of opposite sublattices and A is a constant. The sum over r ensures that the 
spin excitations have net zero momentum as the incoming photon has essentially zero momentum. For 
calculation of the Raman intensity one needs to consider the correlation function of the spin pair operator, 

P, = £S(r).S(r + $) (4) 

r 

Now for the B\ g symmetry we have, 

E sc ~ x + y. 

So that the Hamiltonian assumes the form, 

H R ~ ]T S(r).{S(r+p£)-S(r + ra)} (5) 

r,p=±l 

Now let us consider the diagram representing (Fig. 1) the interaction between the magnons in real space. 
The interaction between the magnons excited at sites ri and ri + fi can be written as, 

Vint = -J]T 5-(n)5 + (n)5-(n + v)S + (n + „) (6) 
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where \x connects to NN sites of the opposite sublattice. The form of the interaction indicates that the 
propagators on either sides of the interaction line are coupled as \x can connect both in the x and y direction. 
But they can be decoupled if one rewrites the interaction term as (see Fig. 1), 



ri L 

+5o 8 (r 1 ),S+(r 1 ) + So 4 (r 1 )S+(r 1 ) 
where the 4 operators are defined in the following way, 

( r ) = \ + -S ± (r + P y)} 



(7) 



S o 2 W = ^{^(r+px) + S ± (r + py)} 
p 

s o^) = lJ2{p s±{ - r+ p^-p s±( - r+ py» 



S o» - \Y J {pS ± {^+vi)+pS ± {v+py)} 

v 

where the +sign and the —sign in the superscript refers to raising and lowering operators respectively. The 
first two operators correspond to B\ g and A\ g symmetry factors respectively. Hence explicitly we write them 
as Sr and for our convenience. In this new notation the interaction term looks like, 

V i ^ = -jY l S-{r 1 )S + {T 1 )\SsJr 1 )S+Jr 1 ) + S^JrJS+Jn) (8) 

ri 

+5 03 (n)5+ (r 1 ) + 5o 4 (r 1 )S+ (n) 

Having written the interaction term we proceed to define the spin-wave propagators in different scattering 
geometries. In Fig. 1 the spin raising and lowering operators at r and r + 8 are denoted by S + (r) and 
Sg (r) and those at at r' and r' + 8' are denoted by S~(r') and Sg (r'). This immediately suggests that 
8 and 8' which connect NN pairs are fixed in space and we can write the spin- wave propagators as, 



X 



-(r,r') = (^G|S-(r)S + (r')|^G> 



(9) 
(10) 



corresponding to the upper line and, 

XB -(r,r>) = (ip G \S B -Jr)S B Jr>)\iP G ) 

corresponding to the lower one in Fig. 1. Now the time-ordered two-magnon propagator in the B\ g geometry 
can be defined as, 

G Bl » = -if ^( t - i, )(fc|rEs-(r,i)S+ s (r I f)S+(r',(')S Bi! (r',t')]|fc) (11) 

J r,r' 
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For calculating the Raman intensity first consider the non-interacting limit of the two-magnon propagator, 



r,r ; 

The configuration averaged two-magnon Raman propagator is given by, 



(12) 



(13) 



We can decouple the two propagators if we neglect the corrections arising due to exchange between the 
magnons. This allows us to write, 

^>) = * / ^E^ + ( r ' r ''^) X^~(r,r'^-^) (14) 

•* r,r' 

Having done the configuration averaging of G°(w) one can proceed to calculate the two-magnon Raman 
intensity which is obtained from imaginary part of G(u>) where, 



G{lo) 



1 + JG H 



(15) 



We can neglect the vertex corrections since they appear only at the second order level and hence proportional 



to the fourth power of the disorder strength i.e. (— ) (in reference [5] this result is shown explicitly for the 
on-site (diagonal) disorder case and calculation of the vertex correction is discussed). So for weak disorder, 
the correction will be insignificant. This enables to write, 



1 + JGoH 



(16) 



Now it is shown in Appendix A that the B\ g symmetry factor viz. (cosg x — cosgj,) 2 corresponding to the 
pure case is obtained by Fourier transforming the non-interacting two-magnon propagator G° Bi (Eq. 12). 
This symmetry factor is important for calculating the two-magnon Raman scattering intensity. It can also 
be noted that due to the form of the interaction i.e. presence of both B\ g and A\ g terms in Vm, off-diagonal 
combinations of different scattering geometries in G°(uj) are possible. In that case the symmetry factor will 
contain contributions from different geometries. We can see it clearly if one considers, 

(/Sjf (Bi g )S~, (Ai g )) which when Fourier transformed yields, 



E 



(cosq x - cos g^) (cos ^ + cosq y )x + (q,u) 



(17) 



But this symmetry term can be shown to be identically equal to zero for all frequencies because of the 
symmetry present in the magnetic Brillouin Zone. Hence the only terms which survive are the ones that are 
diagonal in scattering geometries M . 



III. MODE ASSIGNMENTS 



Here we want to establish a one-to-one correspondence between the pure system and the one in presence 
of hopping disorder. For doing this we calculate the magnon modes using the exact eigenstate method 
for a 10 x 10 system in absence of disorder in the strong-coupling limit (U/t = 200) . This procedure of 
obtaining the magnon modes numerically has been discussed in reference [4]. These modes thus obtained 
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are assigned a pair of mode numbers n x and n y which are related the wavcvectors q x and q y via the relation 
q x / y — ^-n x / y in the translationally invariant (pure) case. The numbers n x and n y are chosen such that 
the corresponding q x and q y are restricted to the upper-half (first two quadrants) of the Brillouin Zone (BZ) 
where q x is allowed to have both positive and negative values (ranging from — ir to +ir) whereas q y is assigned 
only positive values (from to ir) . The other half of the BZ does not need to be considered as the q x and 
q y values (or equivalently the n x and n y ) lying within the region specified (see Fig. 2) are sufficient to label 
all the modes. The Goldstone mode or the lowest energy mode is assigned with the mode numbers n x — 0, 
n y = 0, the second collective mode with n x = 1, n y — (or equivalently n x ~ 0, n y — 1), the third one with 
n x = 1, n v = 1 and so on. Having done with the assigning scheme it is now possible to form the symmetry 
factor for the B lg scattering geometry which is 0^ = [cos^)^ — cos(^-)n y } 2 . This symmetry factor is 
necessary for calculating the two-magnon Raman intensity. The consistency of the mode assigning scheme 
is checked by re-evaluating the magnon energies using the mode numbers via the relation uj n — 2 Jyl — 
where 7„ = ^[cos^)?"^ + cos(^)n y ] and J = ^ = 0.02 (i.e. U/t = 200). The energies thus obtained 
show excellent agreement with the ones calculated from the exact eigenstate method for the pure case. The 
same mode assignments are expected to be carried over for the intermediate coupling regime. We show 
the comparison between the energies calculated numerically from the exact eigenstate method and the ones 
evaluated by using the analytical expression relevant to the intermediate-C/ regime (U jt = 10) viz. 



UJ n = 2J 



e r /2tt\ /2tt\ i 1/2 



(!-7j - ^{6 + 3cos I — I n x cos I — I n y - 9 7 „} 



(18) 



in Table I. We establish the correspondence between the pure and the disordered system by assuming the 
same assigning scheme continues to work for the disordered system. We show the mode energies for the 
disordered system (a — 0.1) for a particular configuration in Table II. We also notice that the degeneracy of 
the modes is lifted due to the presence of disorder and the low energy modes are not significantly affected 
whereas there is a considerable effect towards the higher energy side. To substantiate this result regarding 
the nature of the high energy modes, we have examined the magnon wavefunction for the high-energy modes 
and find that these modes are strongly localised in certain regions of the lattice fl. In reference [3] it is 
also mentioned that the localization of the high-energy magnon modes occurs in those regions of the lattice 
where the locally-averaged hopping strength is significantly higher than the bulk average, so that the local 
magnetization, which goes as ~ U/t is low. This is shown in Fig. 3. 

IV. CALCULATION OF THE TWO-MAGNON RAMAN INTENSITY 

As before, we can express the propagator in terms of it's non-interacting components Gq{uj). The 
mode assignments allows us to write Gq(uj) in a 2 x 2 matrix form with the mode energies and the symmetry 



factors are labeled by the mode numbers n(= y n x + n y )- Hence we can sum up the perturbation series in 
powers of the interaction between the magnons and which yields, 

> G "Mi = Tr|b (19) 

As in reference [7] the above equation is written in terms of the matrix elements, 

[G 0n (oj)} A A = An, 

[G 0n (u)]BB = B n and 

[Go n (ui)]AB = [Go„(w)]sA = C n 



Trn _ A n + B n +2C n + 2J(A n B n -C%) 

Gn[U}) - (l + JA n )(l + JB n )-J*CZ (20) 
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where A ni B n and C n arc defined as 

2 



An — 



n 



2. J 



2 J 



2 J 



T7 ^ ' 



2J 



2 V2J7 
1 /u\ s ' 



1 



1 



V2J/ V 2J / 2 V2J/ 



w — 2uj n ui + 2ui n 
1 1 



ui — 2ui n ui + 2ui n 



where, 



m= 1 - 2 
6= 1 



A 7 ' 
11 t 2 
2 



= 1 



1 



1 



w — 2u;„ to + 2w r , 



{3 + § cos^n*) cos(^n,,) + 7;?} and 



In Fig. 4 we have shown the Raman scattering spectra in the B\ g symmetry as a function of transferred 
photon frequency u> (in units of J) for a system size 10 x 10 in presence of disorder (<r = 0.1) and U/t = 10. 
It is noticed that the linewidth to peak position ratio is even higher than unity, whereas the experimentally 
value is closed to 0.38 |1. This shows clearly that the strength of disorder a (or equivalently the mean square 

deviation {[p^j-) 2 ^) does not need to be as large as 0.1 in order to get an agreement with the experimentally 
observed data for the Raman spectrum. Fig. 5 shows the same plot for more moderate values of a viz 
a = 0.01 and 0.03. The spectral lineshape for the pure case (a = 0) is included for comparison. This plot 
is reminiscent of the magnon DOS shown in Fig. 2 in reference [3]. There we have considered much larger 
strengths of disorder (a — 0.1 and 0.3). It can be observed that the linewidth to peak position ratio increases 
from 0.283 corresponding to the pure case to 0.37 for a = 0.01 and 0.410 for a = 0.03 whereas the position of 
the peak roughly stays the same in all the cases. So it looks like a = 0.01 roughly corresponds to the value 
that is needed for parameterizing disorder strength in order to obtain agreement with the experimentally 
obtained Raman lineshape. Also it is seen that in presence of disorder there is a pronounced asymmetry 
in the higher frequency region of the spectrum and a considerable intensity is observed all the way upto 
6J and even more. This can be explained by arguing that the high energy magnon modes are strongly 
affected due to a cooperative effect arising from local correlations of hopping disorder, which resulted in an 
appreciable one-magnon DOS beyond the maximum energy 2 J for the pure system ||. So it seems that 
a much lesser value of the disorder strength can actually account for the puzzling features of the Raman 
lineshape. Finally, in Fig. 6 we plot the scattering spectrum in the strong coupling limit (U jt = 100) for 
various disorder strengths viz. a = 0.01, 0.03 and 0.1 in order to draw comparisons with the results obtained 
by Nori et al. [jjj since their work deals with the strong coupling limit i.e. U/t —> oo. The linewidth to peak 
position ratio in this case comes out to be 0.333, 0.406 and 1 corresponding to a — 0.01, 0.03 and 0.1. Also 
there is a pronounced asymmetry with respect to the two-magnon peak in the higher frequency regime. To 
obtain an idea about the asymmetry we measure the quantity I{u> = AJ) / I max (uj) . This has a value close to 
0.4 which agrees very well with the corresponding quantity obtained from the experiments done on La2Cu04 
[g. Nori et al. obtained agreement with experiments for <jj ~ 0.5 J (they have defined crj in units of J) 
which is much more than what we have considered. This can be explained by noting that their Quantum 
Monte Carlo (QMC) calculations of the Raman spectrum for a S = 1/2 Heisenberg antiferromagnet in 
the Big scattering geometry reveals a much lesser linewidth to peak-position ratio (~ 0.176) since in their 
calculation the intrinsic magnon damping Q| has not been taken into account (remember we obtained 0.283 
for the corresponding quantity |7j which is much closer to the experimental value 0.38 than that of Nori et 
al). This explains the reason why they need a much bigger width (<r) in the distribution of the exchange 
constant J so as to have an agreement with experiments. In fact, the evaluation of the different spectral 
moments of the Raman lineshape due to Singh et al. || shows the ratio of the second moment to that of 
the first one (which is like the linewidth to peak-position ratio) as 0.23 which tallies very well with the 
experimental value 0.28. 

In conclusion, we have studied the effects of hopping disorder on the two-magnon Raman intensity using the 
magnon energies calculated numerically using the exact eigenstate analysis. The energies for the disordered 
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case are assumed to have a one-to-one correspondence with the ones in the pure case. This correspondence is 
established by assigning mode numbers n x and n y which are related to the momentum labels via the relation 
Qx/y — {^j;) n x/y The validity of the assigning scheme is checked for the pure case by re-evaluating the mode 
energies using the mode numbers and hence comparing them with the ones obtained numerically. Once 
the consistency of the assigning scheme is checked, the numbers can be used to determine the symmetry 
factors which are required for calculating the Raman intensity in different scattering geometries. For the 
Raman lineshape in the B\ g in presence of a small amount of disorder (<r as small as 0.01) symmetry it 
is observed that the lineshape broadens and alongwith there appears a substantial asymmetry towards the 
higher frequency side while the low-energy modes are weakly affected. It is also observed that the spectral 
intensity persists even beyond 6 J where the joint magnon DOS diverges at 4J. We also obtained a good 
agreement for the linewidth and asymmetry of the two-magnon profile observed experimentally from the 
Raman studies on cuprate insulators for a = 0.01. This clearly reflects the role of zero-point fluctuations of 
the lattice to explain the anomalous features appearing in Raman spectrum for the cuprate antiferromagnets 
via a highly asymmetric magnon energy renormalization arising due to correlation of hopping disorder. 



APPENDIX A 



Here we find an expression for the symmetry factor appropriate to the B\ g scattering geometry. The 
Flcury-Loudon Hamiltonian for the B\ g scattering geometry is written as (see Eq. 5), 

H = A S(r).{S(r+px)-S(r+py)} (21) 

r,p— ±1 

The two-magnon propagator in the non-interacting limit for this case is, 

GoH=* J ^'£x- + (r,T',Sl 1 )xt-(T,?,w-w 1 ) (22) 

r.r' 

where r, r' denote lattice sites and f = r + 5, r' = r' + 5', where 5/5' connect to the sites neighbouring to r 
and r' in all the four directions i.e. 6/5' can be px and py with p = ±1. Also the spin raising and lowering 
operators forming the propagator is shown in the Fig. 7. 

Clearly consists of 16 terms with 5 and 5' having all possible combinations which are enumerated 
below (the energy dependence is dropped temporarily). 
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-x + - 


(r 


+ y, r ' 


— x) 
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To obtain the symmetry factor one has to Fourier transform all the 16 terms in Gq{u). As an example, the 
calculations for the first two terms are shown below, 

(1) The first term is, 
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i f ~^^2 X + ( r ' r ')x + (r + x,r' + x) 

r,r' 

Taking the Fourier Transform, 

* 1 ^EE^W ei! - (r " r,) ^(«> i? '- (r " r,) 

r,r' 

Which can be written as, 



(2) And the second term, 



1 f ~^J2 X + ( r ' r ')x + (r + x,r'-x) 



Fourier Transforming the above quantity yields, 

1 J ^EE^ + W e4 "- (M ' ) x + "(9 , )^'- (r - r '- 2£) 



r,r' q,q r 



Which is equivalent to, 



* j '^J2 x + (^ x+ (-i) 



e 2ig£ 



Likewise we can calculate the contributions from the other 14 terms and all of them combine to give, 

-7j^-^2x~ + (q,^i)x + ~(-q,u - ^[(cosfe - cosg y ) 2 ] (24) 

Hence the symmetry factor used in the B\ g geometry is (cosq x — cos,q y ) 2 . By taking the appropriate signs 
(positive) for the spin operators corresponding to the lower line in Fig. 7 for the A\ g geometry, it is trivial 
to show that the symmetry factor is {cosq x + cosq y ) 2 . 
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S(r) 



+ 



S(n) S(n) 



S(r') 




r 




r' 



r+ 5 



r' + 5' 






FIG. 1. Interaction between two magnons is shown where the bold lines represent the two magnons propagating 
from neighbouring sites r and r + 8 to r' and r' + 8' . The wavy lines denotes the interaction between the magnons at 
sites ri and ri + fi where fi connects to the nearest neighbouring sites. The creation and destruction of magnons at 
these points are shown by S + and S~ where the subscripts B\ g and Ai g correspond to different scattering geometries 
defined in the text. 
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FIG. 2. The region in the Brillouin Zone (BZ) is shown where the dots denote the mode numbers n x and n y 
chosen so as to assign labels to to the mode energies calculated. A few of the mode numbers are indicated. The first 
two quadrants are required for assigning numbers to all the modes. 
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FIG. 3. The amplitude for the highest energy magnon mode (below) is shown for a 10 x 10 lattice for U/t = 10 in 
presence of disorder a — 0.1. Also we have shown the magnetization profile (top) which shows a dip at places where 
the locally averaged hopping strength is maximum. 






1 2 3 4 5 6 



©/J 

FIG. 6. The two-magnon Raman lineshape is shown for the strong coupling case (U /t — 100) in order to compare 
with the results of Reference [1]. The values taken for disorder are a = 0.01, 0.03 and 0.1. 
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+ 




r+5 r ' + 5 



^ <r 



"2 (Sr +x- Sr + y) "^(Sr'+x S r '+y) 

FIG. 7. The two-magnon Raman propagator in the non-interacting limit corresponding to the Bi g symmetry is 
shown. 
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TABLE I. The mode energies are obtained using the expression relevant to the intermediate-^ regime (U/t = 10 
or A/t = 5) for a 10 x 10 system where n x / y = {^;)<lx/ y are the mode numbers specified in the leftmost column. Also 
we have shown the energies calculated numerically using the exact eigenstate method for comparison. The symmetry 
factors for both Bi g and A\ g scattering geometries are obtained using the mode numbers. 



Mode Nos. (n x , 


n y ) 


ui n (numerical) 


u n (formula) 






(-5,0) 




0.7548 


0.7505 


4.000 


0.000 


(4,0) (0,T 


(-4,0) 


(-5,1) 


0.7438 


0.7408 


2.618 


0.000 


(4,1) (1,4) 


(-4,1) 


(-1,4) 


0.7396 


0.7356 


0.382 


0.000 


(3,2) (2,3; 


(-3,2) 


(-2,3) 


0.7150 


0.6931 


3.272 


0.037 


(3,1) (1,3; 


(-3,1) 


(-1,3) 


0.7030 


0.6927 


1.250 


0.250 


(4,2) (2,4; 


(-4,2) 


(-2,4) 










(3,0) (0,3; 


(-3,0) 


(-5,2) 


0.6882 


0.6793 


0.000 


0.382 


(2,2) (3,3; 


(-2,2) 


(-3,3) 


0.6765 


0.6696 


1.713 


0.478 


(2,1) (1,2; 


(-2,1) 


(-1,2) 


0.5945 


0.5826 


0.250 


1.250 


(4,3) (3,4; 


(-4,3) 


(-3,4) 










(2,0) (0,2; 


(-2,0) 


(-5,3) 


0.5469 


0.5458 


0.477 


1.714 


(1.1) (4,4; 


(-1,1) 


(-4,4) 


0.4199 


0.4102 


0.000 


2.618 


(1,0) (0,1; 


(-1,0) 


(-5,4) 


0.3061 


0.2987 


0.036 


3.273 




0,0) 




0.0000 


0.0000 


0.000 


4.000 





TABLE II. Here the assignment of the mode numbers n x and n y is shown for the mode energies calculated 
numerically for both pure and disordered systems (a — 0.1) of size 10 x 10 and for U/t = 10. 



Mode Nos. (n x , n y ) 


io n (pure case) 


u>„ (disordered case) 


(-5,0) 


0.7548 


1.3006 


(4,0) (0,4) (-4,0) (-5,1) 


0.7438 


1.2222, 1.2024, 1.1960, 1.1346 


(4,1) (1,4) (-4,1) (-1,4) 


0.7396 


1.0578, 1.0178, 1.0137, 0.9955 


(3,2) (2,3) (-3,2) (-2,3) 


0.7150 


0.8557, 0.8414, 0.8291, 0.7934 


(3.1) (1,3) (-3,1) (-1,3) 

(4.2) (2,4) (-4,2) (-2,4) 


0.7030 


0.7854, 0.7762, 0.7698, 0.7641 
0.7228, 0.7018, 0.6935, 0.6624 


(3,0) (0,3) (-3,0) (-5,2) 


0.6882 


0.6614, 0.6554, 0.6364, 0.6197 


(2,2) (3,3) (-2,2) (-3,3) 


0.6765 


0.6020, 0.5773, 0.5668, 0.5411 


(2,1) (1,2) (-2,1) (-1,2) 


0.5945 


0.5251, 0.5167, 0.4926, 0.4784 
0.4527, 0.4361, 0.4256, 0.4151 


(2,0) (0,2) (-2,0) (-5,3) 


0.5469 


0.3969, 0.3711, 0.3595, 0.3424 


(1,1) (4,4) (-1,1) (-4,4) 


0.4199 


0.3320, 0.3170, 0.2909, 0.2778 


(1,0) (0,1) (-1,0) (-5,4) 


0.3061 


0.2390, 0.2118, 0.1801, 0.1643 


(0,0) 


0.0000 


0.0000 
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